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In this paper we study nuclear effects in the neutrinoproduction of pions. We found that in a 
Bjorken kinematics, for moderate xb accessible in ongoing and forthcoming neutrino experiments, 
the cross-section is dominated by the incoherent contribution; the coherent contribution becomes 
visible only for small |t| < 1/ R\, which requires xb <0.1. Our results could be relevant to the kine¬ 
matics of the ongoing MINERvA experiment in the middle-energy (ME) regime. We provide a code 
which could be used for the evaluation of the i/DVMP observables using different parametrizations 
of GPDs and different models of nuclear structure. 


I. INTRODUCTION 

Today one of the key objects used to parametrize the nonperturbative structure of the target are the generalized 
parton distributions (GPDs). For kinematics where the collinear factorization is applicable Ell , they allow to evaluate 
cross-sections for a wide class of processes. Right now all the information on GPDs comes from electron-proton and 
positron-proton measurements done at JLAB and HERA, in particular the deeply virtual Compton scattering (DVCS) 
and deeply virtual meson production (DVMP) 0MEM3- A planned CLAS12 upgrade at JLAB m and ongoing 
experiments at COMPASS [Ts| will help to improve our understanding of the GPDs, and in particular the ability to 
polarize both the beam and the target will allow to measure a large number of polarization asymmetries, providing 
various constraints for phenomenological GPD parametrizations. However, in practice the extraction of GPDs from 
modern experimental data is still aggravated by uncertainties, such as large BFKL-type logarithms in next-to-leading 
order (NLO) corrections [19] at HERA kinematics, higher-twist components of GPDs and pion distribution amplitudes 
(DAs) at JLAB kinematics [2flfj23| . and vector meson DAs in the case of p- and <fi -meson production. 

From this point of view, consistency checks of the GPD extraction from experimental data, especially of their flavor 
structure, are important. Earlier we proposed to study the GPDs in deeply virtual neutrinoproduction of pseudo- 
Goldstone mesons (w, K, rf) [24| with high-intensity NuMI beam at Fermilab, which recently switched to the so-called 
middle-energy (ME) regime [ 23 , with an average neutrino energy of about 6 GeV. The i^DVMP measurements with 
neutrino and antineutrino beams in this kinematics are complementary to the electromagnetic DVMP measured at 
JLAB. In the axial channel, due to chiral symmetry breaking, we have an octet of pseudo-Goldstone bosons, which act 
as a natural probe of the flavor content. Due to the V — A structure of the charged current, in i/DVMP one can access 
simultaneously the unpolarized GPDs, H , E, and the helicity flip GPDs, H and E. Besides, using chiral symmetry 
and assuming closeness of pion and kaon parameters, the full flavor structure of the GPDs may be extracted. We 
found [2fil | that the higher-twist corrections in neutrino production are much smaller than in the electroproduction, 
which gives an additional appeal to the neutrinoproduction channel. 

Unfortunately, in modern neutrino experiments, for various technical reasons, nuclear targets are much more fre¬ 
quently used than liquid hydrogen. By analogy with neutrino-induced deep inelastic scattering (izDIS) on nuclei, 
one can expect that izDVMP on nuclei could be sensitive to many nuclear phenomena such as shadowing, antishad¬ 
owing, EMC-effect and Fermi motion. In inclusive processes, all these effects give contributions of order <10% in 
the 0.1 < x r 5. 0.8 region relevant for the ongoing and forthcoming i/DVMP experiments [27]. However, there are 
indications [281] that in the off-forward kinematics (t ^ t m i n ) they could be enhanced. For this reason, in order 
to be able to test reliably various GPD models, one should take into account nuclear effects. We study them in a 
handbag approach, since in the regime of xr > 0.1, relevant for the current and forthcoming ^DVMP experiments, 
multiparticle corrections should be negligible. 

The paper is organized as follows. In Section [TT] we discuss the framework used for evaluation of nuclear effects. In 
Section HTTl for sake of completeness we list briefly the parametrizations of GPDs used for our analysis. In Section HVl 
we present numerical results and draw conclusions. 


II. NUCLEAR EFFECTS 

There are two types of processes on nuclear targets, coherent (without nuclear breakup) and incoherent (with 
nucleus breakup into fragments). The former contribution is enhanced due to coherence as ~ A 2 , but this effect is 
relevant only at very small values of t: At larger values of \t,\ ~ 1/r^, where ta is the nuclear radius, this contribution 
vanishes rapidly and eventually gets covered by the incoherent contribution. 
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The typical values of xb accessible in the modern and forthcoming ^DVMP measurements are xb > 0.1, and for 
this reason one can neglect multinucleon coherence effects and describe the process by single-nucleon interactions. 
Combining this with the weak binding of the nucleons inside nuclei, we may use the Impulse Approximation (IA) and 
write the amplitude of the process as, 

(E 

(!) 

where k — A/2 and k + A/2 are the momenta of the incoming and outgoing nucleons respectively, p p and p n are 
the density matrices of the protons and neutrons inside nuclei, and A ppl are the amplitudes of the process on free 
protons and neutrons [29^ . In (JT]) we ignore a poorly known and essentially model-dependent contributions of the 
so-called non-nucleonic degrees of freedom, which are sometimes added to the rhs of ©• Also, we don’t include the 
contribution of processes in which a final nucleus remains in an excited isomer A* state: we expect that such processes 
are suppressed both at large-t (due to the nuclear formfactor) and small-t (due to additional factor ~ t n in multipole 
transitions between different shells). 

In a Bjorken kinematics region, a collincar factorization theorem tells us that the scattering amplitude both on 
the nucleons and nuclei has a form of the convolution of the GPD of the baryon H A with a process-dependent hard 
coefficient function C{x,£) 1 , 


Aoh = / d 3 k 



r A r A 


d 3 k p■ 



A:oh ~ J dxC (x, £) Ha (x, £, t ), 

which, combined with yields a convolution relation for the GPDs of the nucleus 2 , 

H q/A (x, £,t)= [ —H p/A (y, f, t) H p {-, ~,t) + [ —p n (y, £, t) H n/A -, t) , 
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where y is the light-cone fraction of the nuclear momentum carried by the nucleon, and we introduced the so-called 
light cone nucleon distributions H p / A , H u /a related to the densities p Pi „ as 


f 2 A ± Aj_ 

Hi/ A {y, t) = mjv / d k±pi [ m N {y + 0, fc-L -—; m N (y - ^), k± + — 


, i=p,n. 


(4) 


The two equations m may be schematically illustrated with the so-called double handbag diagram in the left 
pane of the Figure [U as a two-stage process. This approximation has been used e.g. in |28 h 33| . and describes eA 
data reasonably well. 

For the incoherent processes, we may assume completeness of the final states (the so-called closure approximation), 
and using unitarity, as schematically shown by the diagram in the right pane of the Figure ©, get a similar expression 
for the cross-section of the process ;2!)|, 

CTincoh = J d 3 k ^2 pi (k, kj ai (k, k + A, q'j & J ^ H i/ A(y,0,0)a l (jt, k + A, q^j (5) 

i—p,n ^ i=p,n ' ^ A 

where P A is the momentum of the nucleus, and the last equality in ([5J is valid in a collinear approximation. 

Since the binding energy of a nucleon in the nucleus is very small compared to a mass of the free nucleon, the 
distributions H p / A , H n / A are strongly peaked functions, and in the first approximation may be approximated as |28l | 
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1 Explicit expressions for the leading twist coefficient functions for various ^DVMP processes may be found in [2^1 

2 In what follows we assume for the sake of simplicity that the spin of the nucleus is zero, which is true for most frequently used nuclear 

targets like 12 C, 40 Ca, 40 Ar, 56 Fe, 132 Xe. 
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Figure 1: (color online) Left: Double handbag diagram for the amplitude of the coherent pion production on the nucleus. 
Right: Closure approximation and its relation to distribution of the cross-section of the incoherent process. 


where Z is the atomic number, A is the mass number, and the parameter a « ~ 200 MeV is the Fermi 

momentum inside the nucleus, which controls the width of the distribution. In the extreme limit a —>• 0, the product 
of the exponent in (TO and a prefactor yja/ir reduce to a (5-function, and instead of a convolution we end up 
with a mere sum of amplitudes (for the coherent case) or cross-sections (for the incoherent case) on separate nucleons. 
However, such a factorized form is an oversimplification, since it cannot describe the A-dependence of the first moment 
of the so-called D -term ^(0), for which there are estimates based on very general assumptions |34j. 

A more realistic approach is to use the functions H p /a, H n /A, evaluated in the shell model of the nuclear structure. 
One of the most popular choices for the evaluation of the nucleon dynamics inside a nucleus is a QHD-I model proposed 
in (35| - |37| . The lagrangian of this model, in its simplest form describes the interaction of the nucleons with effective 
vector and scalar fields, 

C = $(id-M- g v V + g s ^j (d^d^) - ±V liV V> u/ + ^-V^. (8) 

where we used a shorthand notation V pv = and <j> are the fields of vector and scalar mesons respectively. The 

mean field models based on a Lagrangian of type (JHJ), have been successful in the description of various characteristics 
of nuclei. The simplest version of the model used in this work consists of baryons and isoscalar scalar and vector 
mesons. The pseudoscalar pion degrees of freedom are neglected because their contribution to the ground state of 
0 + -nuclei essentially averages to zero !3(i|. In the literature one may find extensions of the model Q, which have 
additional mesonic degrees of freedom and give better quantitative description of nuclei, especially with nonzero spin 
and isospin. The corresponding explicit expression for the distribution functions H p /a, H u /a were calculated in the 
model (0 in [28], yielding 


H Pi n/A {y 5 0 



^, fc_L + 


l+'&i 




(9) 


where <!>i is the wave function of the nucleon inside the nucleus, the summation index i runs over the proton or neutron 
shells respectively. 


III. GPD PARAMETRIZATION 

For numerical estimates of the nuclear effects, one should use a particular parametrization of GPDs available 
from the literature ?|. fl3l [38j - l44| . For the sake of definiteness^ in what follows we use the parametrization of Kroll- 
Goloskokov (H, 46T which succeeded to describe HERA [47] and JLAB [3Sj, 05|,[46] data on electroproduction of 

different mesons, and therefore it should provide a reasonable description of neutrino-induced DVMP. The parametriza¬ 
tion is based on the Radyushkin’s double distribution ansatz, in which the skewness is introduced separately for sea 
and valence quarks, 


— H-va.1 (*G t) T H sea (x, £, t) 


(10) 
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Figure 2: (color online) Ratio (1131) for coherent (left) and incoherent (right) 7r -production for several nuclei. In the lower 
parts of each figure we show the asymmetry m, which measures the difference between 7r + and n production. As explained 
in the text, for the first three nuclei it is exactly zero, so for the sake of legibility we don’t show those curves. 
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and q va i and q sea are the ordinary valence and sea components of PDFs. The coefficients bi, oci, as well as the 
parametrization of the input PDFs q(x), A q(x) and pseudo-PDFs e(x), e(x) (which correspond to the forward limit 
of the GPDs E, E ) are discussed in [ 4 ^, H^|. The unpolarized PDFs q(x) are adjusted to reproduce the CTEQ 
PDFs in the limited range 4 < Q 2 < 40 GeV 2 . The j/DVMP cross-sections on free protons and neutrons have been 
evaluated with this parametrization in our previous papers [24j, [26| . 


IV. RESULTS AND DISCUSSION 


In order to quantify the size of the nuclear effects, we consider a ratio 


Ra 


da a I dt du dQ 2 

Z da p /dt dv dQ 2 + (A — Z) da n /dt dv dQ 2 


(13) 


which takes into account differences in isotopic content of different nuclei. For the case of self-conjugate nuclei, 
this ratio up to a coefficient A/2 coincides with a deuteron-normalized cross-section 3 used in the presentation of 
experimental data. 


The nuclear effects in the deuteron are small. The shadowing corrections are also negligible since we consider the kinematics xg > 0.1. 
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Figure 3: (color online) f-dependence of the ratio 11131) for coherent (left) and incoherent (right) 7r -production for several nuclei. 


In the left pane of Figure (|2]l. we have plotted the ratio Ra da for several spin-0 nuclei. In the regime of small-a ;b 
the ratio Ra is close to A due to the coherence of contributions of separate nucleons. For higher values of xb : the 
behavior of the cross-section is similar to that of a formfactor: it decreases rapidly and has nodes, with average 
distance between the nodes ~ 1 /ta, where ta is the nuclear radius. However, the positions of the nodes do not 
coincide with those of a nuclear formfactor, a type of behavior which cannot be reproduced by a simple model (1517b . 
In neutrino experiments this kinematic region is hardly accessible experimentally, since it is covered by the incoherent 
contribution if the final nucleus breakup is not detected (see the right pane of the same Figure). In the lower pane of 
each figure, we’ve shown the asymmetry 
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— 7 ) - da, 
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7r — n) ^ 


(14) 


which is sensitive to an isospin-1 GPD combination H u — H d . For self-conjugate nuclei, this asymmetry is exactly 
zero since in the model [36j, |37| the difference between proton and neutron distributions is negligible. For 90 Zr and 
208 Pb, the asymmetry 1151) in general is small and does not exceed 10%, although increases slightly near the nodes of 
7r + and 7r~ 4 . 

The f-dependence of the cross-section is shown in Figure [3] At small-f, the coherent ratio Ra scales as Ra oc 
A exp(t mirl (j:s)r^/6), where ta is the nuclear radius, but decreases rapidly at large-t. The incoherent cross-section 
shown schematically in the right pane of Figure [51 is close to unity, as expected. Its suppression at small-|t| ~ \t m i n \ 
comes from < 1 and onshellness conditions in the convolution integral in ([3]). 

In order to understand the sensitivity of the cross-sections to a choice of GPD parametrization, in the left pane of the 
Figure [4] we compared the predictions of Kroll-Goloskokov model discussed in Section ITTT1 with a simple zero-skewness 
model H q = q(x)F]y(t). As we can see, there is an up to a factor of two difference between the two models. 

Frequently the targets in neutrino experiments are organic scintillators with a general atomic structure CH n . As 
one can see from the right pane of the Figure [H in the region xb ^ 0.3 there are two dominant contributions, from 
hydrogen atoms and from incoherent cross-sections, which cannot be separated unless a final nucleus is detected. A 
coherent cross-section is strongly suppressed in this kinematics. 


V. CONCLUSIONS 

In this paper we studied the nuclear effects in the coherent and incoherent pion production. We found that the 
former has a complicated structure, with coherent enhancement in the region of small-a ;b , small-f, and strong nuclear 


4 The nodes of 7r + and n don’t exactly match due to differences in proton and neutron distributions p p , p n . 
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Figure 4: (color online) Left: Comparison of coherent cross-sections evaluated in Kroll-Goloskokov (KG) and Zero Skewness 
(ZS) models. Right: Comparison of different contributions to pion production on atomic targets CH„. 


suppression outside this kinematics. Similar to a formfactor, the leading twist contribution has nodes. For the 
incoherent case, the nuclear dependence is quite mild, and for |t| > 3\t m i n (xB)\ the nuclear effects are negligible, i.e 
the full cross-section is a mere sum of contributions of separate nucleons. This is a model-independent result, and 
from a practical point of view, this allows to get rid of extra uncertainties related to nuclear structure. Our approach 
is applicable in the regime xb > 10~ 1 2 * 4 5 6 7 8 9 10 11 12 13 . For smaller values of Xb, this picture is modified due to coherence [48| and 
saturation [49] effects. For further practical applications, we provide a code, which can be used for the evaluation 
of nuclear cross-sections with different parametrizations of GPDs and models of nuclear structure. The modular 
structure of the code allows to easily consider different parametrizations of GPDs and nuclear distributions. For 
illustration, we provide with this package the libraries for the Kroll-Goloskokov GPD model and the QHD-I nuclear 
structure model distribution used in this paper. Also, we provide detailed instructions how to build and use new 
libraries. 

Finally, we would like to stop briefly on the recent results of the MINERvA collaboration [29]. Albeit those results 
are for coherent pion production on nuclei, here we do not make any comparison, because the conditions of applicability 
of collinear factorization ( Q 2 rn 2 N ) are not met in that kinematics. Models based on extrapolation of the Adler 
relation [50] are more appropriate for that kinematics. 
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